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I. INTRODUCTION 

In general Relativity (GR), the spacetime singularity forms due to gravi- 
tational collapse of a massive astrophysical object. According to singularity 
theorems [1] , the singularities exist in the form of either future or past incom- 
plete spacelike geodesies. In other words, trajectories of material particles 
(photons) will come to a sudden end (disappear from the spacetime). The 
singularity theorems fail to provide practical information on the nature of 
the singularity. 

It has been an interesting open problem to know the final fate of the 
gravitational collapse. Penrose [2] suggested that the spacetime singularity 
is always covered by the event horizon. This is known as the Cosmic Cen- 
sorship Hypothesis (CCH), which has no mathematical or theoretical proof. 
The spacetime singularity would be a black hole or a naked singularity, de- 
pending on the initial data and the equation of state. Many efforts have been 
made to prove or disprove the CCH. For this purpose, Virbhadra et al. [3] 
introduced a new theoretical tool using gravitational lensing. Virbhadra [4] 
used gravitational lensing to find an improved form of the CCH. 

Oppenheimer and Snyder [5] are pioneers who investigated dust collapse 
by taking the static Schwarzschild spacetime as an exterior spacetime and 
a Friedmann like solution as an interior spacetime. They concluded that a 
black hole is the final fate of gravitational collapse. Markovic and Shapiro 
[6] generalized this work by taking a positive cosmo logical constant. Lake 
[7] extended it for both positive and negative cosmological constants. Sharif 
and Ahmad [8] extended spherically-symmetric dust gravitational collapse 
with a positive cosmological constant to a perfect fluid. 

Recent developments in string theory and other field theories indicate that 
gravity is a higher-dimensional interaction. It would be worthwhile to study 
gravitational collapse and singularity formation in higher dimensions. Ghosh 
and Banerjee [9] have studied dust collapse with five- dimensional Tolman- 
Bondi spacetime. In a recent paper, Sharif and Ahmad [10] studied the gravi- 
tational collapse of a perfect fluid in 5D for spherically symmetric spacetimes. 

The behavior of the electromagnetic field in a gravitational field has been 
the subject of interest for researchers over the past decades. The inclusion 
of an electromagnetic field in gravitational collapse predicts that the gravi- 
tational attraction is counterbalanced by the Coulomb repulsive force along 
with the pressure gradient [11]. Sharma et al. [12] concluded that the elec- 
tromagnetic field affects the values of the red-shifts, the luminosity, and the 
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masses of relativistic compact objects. Nath et al. [13] demonstrated that 
the electromagnetic field reduces the pressure during gravitational collapse. 
Recently, we [14] found that the electromagnetic field, along with the matter 
field, in gravitational collapse increases the rate of collapse. 

It would be interesting to study the gravitational collapse of a perfect 
fluid in 5D with an electromagnetic field. The junction conditions between 
the spherically-symmetric spacetimes are discussed using the Israel criteria 
modified by Santos. A closed form of the exact solution of the field equations 
in 5D exists for marginally bound (F(r) = 1) and non-marginally (F(r) ^ 1) 
bound cases. This distinguishes the 5D case from the AD, in which only a 
marginally bound solution is possible. The main objectives of this work are 
to study the effects of an electromagnetic field on the rate of collapse in 5D 
and to see whether or not the CCH is valid in this framework. 

The plan of the paper is as follows: In the next section, the junction 
conditions are given. In Section III, we find the 5D spherically symmetric 
perfect fluid solution of the Einstein field equations with a positive cosmo- 
logical constant in the presence of an electromagnetic field for the marginally 
bound (F(r) = 1) and the non-marginally bound (F(r) ^ 1) cases. The 
apparent horizons and their physical significance are presented in section IV. 
We conclude our discussion in the last section. The geometrized units (i.e., 
the gravitational constant G—l and speed of light in vacuum c =1 so that 
M = and k = = 8ir) are used. All the Latin and Greek indices vary 
from to 4; otherwise, the indices will be given. 



II. JUNCTION CONDITIONS 

A timelike AD hypersurface £ is taken such that it divides a 5D space- 
time into two 5D manifolds, V~ and V + , respectively. The 5D spherically 
symmetric spacetime is taken as an interior manifold V~ [15] 

ds 2 _ = dt 2 - X 2 dr 2 - Y 2 (d9 2 + sin 2 9d(j) 2 + sin 2 9 sin 2 0# 2 ), (1) 

where X = X(t,r) and Y = Y(r,t). For the exterior manifold V + , we take 
the 5D Reissner-Nordstrom de-Sitter spacetime 

ds 2 + = ZdT 2 - \-dR 2 - R 2 (d6 2 + sin 2 9d(j) 2 + sin 2 9 sin 2 0# 2 ), (2) 

Zj 

where 

M Q 2 A , 
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M and A are constants and Q is the charge. The Israel junction conditions 
modified by Santos [16, 17] are 

1. The continuity of the first fundamental form over E gives 

(ds 2 _h = (ds 2 + h = dsl. (4) 

2. The continuity of the second fundamental form over £ gives 

[K ij \=K+-Kr j = 0, (i, j = 0,2,3,4), (5) 
where Kij is the extrinsic curvature defined as 

± ± a 2 24 

correspond to the coordinates on S and x a ± are the coordinates for V^. The 
Christoffel symbols are found from the interior or the exterior spacetime 
and are components of the outward unit normals to £ in the coordinates 
x a ± . 

The equations of the hypersurface £ in terms of the coordinates of the 
interior and the exterior spacetimes are given as [18] 

h-(r,t) = r — r s = 0, (7) 
h+{R,T) = R-R s {T) = 0, (8) 

where is a constant. Using Eqs. (171) and (JED in Eqs.flTJ and (j2D respectively, 
it follows that 

(ds 2 _)x = dt 2 - [Y( rj: ,t)} 2 (d9 2 + sin 2 0# 2 + sin 2 sin 2 0# 2 ), (9) 

- i? s 2 (rf6' 2 + sin 2 ^0 2 + sin 2 0sin 2 ( /)# 2 )- (10) 
For T to be a timelike coordinate, we assume that 



From Eqs. ®, © and JTOj, it follows that 



i?s = ^E,t), (12) 
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The outward unit normals to £ in M and M + are 

n- = (0,X(r s ,t), 0,0,0), (14) 

n+ = (-i? E ,T, 0,0,0). (15) 

The components of the extrinsic curvature Kf- turn out to be 

K 00 = 0, (16) 

YY' 

= esc 2 6*^33 = esc 2 9 esc 2 4>K^ 4 = )sj (17) 

= esc 2 9K+ = esc 2 esc 2 cj)K+ = {ZRT) S , (19) 

where dot and prime indicate differentiations with respect to t and r, respec- 
tively. Continuity of the extrinsic curvature gives 

K+=0, K+=K 22 . (20) 

When we make use of Eqs. (H6])-(|20D along with Eqs. ©, (JT2) and the 
junction conditions turn out to be 

(XY> - XY% = 0, (21) 

M=(Y 2 + ^-^ + Y 2 Y 2 - y 2 ^) s . (22) 



Equations. (1121) . ({TBI and Eqs. (1211 . (1221) are the necessary and the sufficient 
conditions for matching of the interior and the exterior regions. 

III. SOLUTION OF THE EINSTEIN FIELD 

EQUATIONS 

The Einstein field equations with a cosmological constant are 

G$-KS£ = K{Tj! + E%). (23) 
The energy-momentum tensor for a perfect fluid is 

TZ = (p + P )u»u v -p5^ (24) 
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where p is the energy density, p is the pressure, and w M =(l,0,0,0,0) is the 
5- vector co-moving velocity. is the energy-momentum tensor for the elec- 
tromagnetic field given by 

E$ = ^-{-F^F VU + -XF 5uJ F s "). (25) 

First, we solve the Maxwell field equations 

Fpv = 4>u,^ ~ <f>p,u, (26) 
F^„ = -4ttJ^, (27) 

where M is the five potential and J M is the five current. Because we have 
taken the charged fluid in a co-moving coordinate system, the magnetic field 
will be zero in this case. Consequently, the five potential and the five current 
can be taken as 

^ = (0(t,r),O,O,O,O), (28) 
J» = a5l (29) 

where a is the charge density. 

We treat \x and v as local coordinates for the solution of Maxwell field 
equation, Eq. (127]) . Thus, the non-zero components of the field tensor are 
evaluated by using Eqs. ( 1261) and (|28~|) as 

F tr = ~F rt = (30) 



Also, Eqs. ((27D and (J29J) yield 

1 d 2 (j) defy X' 

X dr 2 dr X 2 

1 d 2 (p X d(j) 



X drdt X 2 dr 
The last equation implies that 



-4vraX, (31) 
0. (32) 



(1£) = *<0, (33) 

where K is an integrating function and gives the electromagnetic field con- 
tribution. Equations (l3"Tj) and (I3"3"j) yield 

K'(r) = -AnaX. (34) 
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The field equations in Eq. (f2"3"j) for the interior spacetime are written as 

t ~ xi ' y X Y ^ Y 2 Y 2 X Y Y 2 

= A + K 2 + 8np, (35) 

3Y' 2 Y Y 2 3 
= -^^ + 3( 7 + ^) + ^ = A-8 7 rp + K 2 , (36) 

g 2 Y" X'Y' Y' 2 Y XY Y 2 X 

6 ~ ~X^^Y~YV~2Y 2 ' + 2 ^Y + XY + 2Y 2 + 2X^ 

+ ^ = A-8np-K 2 , (37) 

G\ = G^ = G s e = A-87rp-K 2 , (38) 
Y' X V 

Gl = -3y + 3- y = 0. (39) 

To solve these equations, we first integrate Eq. f l39|) with respect to t so 
that 

Y' 

X = T , (40) 

where F = F(r) represents the energy inside the hypersurface S. Making 
use of Eqs. f )36|) and (pEDl . we get 



r + ( I). + l^! = ^±*!. (41) 

The energy-momentum conservation equation for matter, T^. u = 0, shows 
that the pressure is a function of t only, i.e., p = p(t). Using this value of p 
in the above equation, it follows that 

Y ,Y. 2 1-F 2 A — 8np(t) + K 2 

y + ( y )+ ^^ = 3 ' (42) 

Here, we consider p as a polynomial in t given by [8] 

P(t)=Po(^r q , (43) 

where T is a constant time introduced in the problem by a re-scaling of t for 
physical reasons and po and g are constants. Further, for simplicity, we take 
q = so that p(t) = p . 
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Inserting this value of p(t) in Eq. (142 p and then integrating with respect 
to t, we have 



Y 2 = F 2 - 1 + |£ + (A - 8np + K 2 )^-, (44) 

where m is an arbitrary function of r and is related to the mass of the 
collapsing system. Using Eqs. (jSj) and ( HU1) in Eq. ( )35l) . we obtain 

ml = + Po )Y 3 Y' - (45) 
6 o 

For physical reasons, we assume the following energy conditions [1] are sat- 
isfied: 

P>0, p + p>0, (46) 
P > 0, -p<P< P- (47) 

Integrating Eq. ( )45i) with respect to r, we obtain 

m(r) = y ^(P + Po)^ 3 V" - l£ KK'Y\ (48) 

where m = due to the finite distribution of the mass at the origin (r = 0). 
The mass function m(r) > because m(r) < is not meaningful. Using Eqs. 
(HOD and flSJ) in the junction condition, Eq. (J22|) . it follows that 

M = ^ + 2m + i(ir 2 - 8vrp c )y 4 . (49) 

This mass is located at the origin of the spherical symmetry and produces a 
gravitational field in the exterior region of the sphere. 

The total energy M(r, t) up to a radius r at time t inside the hypersurface 
£ can be evaluated by using the definition of Misner- Sharp mass [18] given 
by 

M(r,t) = ^Y(l+g^Y,Y u ). (50) 
For the interior metric, it takes the form 

M(r,t) = iF(l + y 2 -(^) 2 ). (51) 
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Putting Eqs. (gQD and flS]) in Eq. (JEE]), we obtain 

M(r,t) = ^ + (A + K 2 - 8np )^. (52) 

Now we solve Eq. (jUj) for the following two cases: 

F(r) = 1, F(r) ^ 1. 
1. Solution with F(r) = 1 

First, we discuss the solution with a positive pressure. For A — Sixpo + K 2 > 0, 
the analytic solutions in closed form can be obtained from Eqs. ( 140]) and (1441) 
as follows: 

12m . i , i . 

Fr '* = ("a q —^)*sinh*a(r,t), 53 

A — 87rp + 

12m .! rr 4m' 2KK' . . , . . 

X = (- — )3 { }sinha r,t 



+ ^)-^' ^ (A+ J-«^ ^WV (A + " 2 a" 8m) > 

x cosh a(r, t)] sinh~ a(r, t), (54) 

where 

»M) = ^S[W- f] . (55) 

Here, io( r ) is an arbitrary function of r and is related to the time of the 
formation of the singularity. In the limit (87rp — K 2 ) — > A, the above 
solution corresponds to the 5D Tolman-Bondi solution [19] 

lim Y(r,t) = [8m(t - t) 2 }\ (56) 

(8irp -K 2 )^A 

m'(tn — t) + 2mtL 
lim X(r,t) = — ^ — r-. 57 

(8-po-^ 2 HA [32m 3 (t -t) 2 ]3 
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2. Solution with F(r) ^ 1 

Integrating Eq. QHJ) with the conditions A — 87rp + k 2 > and F(r) ^ 1, it 
follows that 

w n rr 12 ™ 9(F 2 -1) 2 _i . , . . 

A — 87Tpo + -rf (A — 87TJ9 + K j 



3(F 2 - 1) 



A-87rp + K 2i 

Using Eq. (1581) in Eq. (HOI) , we obtain 

1 rr 12m 9(F 2 -1) 2 ,i . , . . 

X = 2F A — 8n Po + K 2 ~ (A — 8rr Po + if 2 ) 2 } ' Smh ^ * } 
3(F 2 -1) n _i r l 12m 9(F 2 -1) 2 

nil o I tT 7 ). ~ 7a o I TP) \ 9 J 



(58) 



A-87rp + K 2J L 2 L A-8vrj9 + ^ 2 (A - 87rp + ^ 2 ) 2 

12m' 24mKK' 36FF'(F 2 - 1) 

X A - 87rpo + if 2 ~ (A - 8vrp + K 2 ) 2 ~ (A - 87rp + ^ 2 ) 2 

36KK'(F 2 - l) 2 . , . . r 12m 

+ (A - 8, Po + K 2 f ] Smh ^ t} + { A - Svrpo + K 2 

9(F 2 -1) 2 i 2KiT 

(A - Sttpo + i^ 2 ) 2 1 2 1 6 (A - Svrpo + K 2 ) I ^ " } 

/ 2(A-8ttpo + ^ 2 ) n i n , , x , 6(W 2 -1)KK' 
+ (- ^ ^)2t '}cosha(r,t) + { 



3 ' u J v ' ' L (A - 87rpo + K 2 ) 2 

§FF' 

}}, (59) 



A - 87rpo + if 2 



where a(r,t) is given by Eq. (155 p . If the dominant energy condition holds 
(i.e., negative pressure), the restriction A — 8%p + K 2 > is removed, and 
a solution is possible for all values of A, p , and K. Equations (158")) and ( |5*0|) 
represent the non-marginally bound solution corresponding to F(r) ^ 1. 
One can easily verify the marginally bound solution given by Eqs. ( 153]) 
and ( 154} by substituting F(r) = 1 into Eqs. ( 158} and ( l59l) . The non- 
marginally bound solution is impossible in a 4D gravitational collapse with 
an electromagnetic field [14]. For K = 0, Eqs. (E3}, ((53}, (EH}, and ((59} 
reduce to the marginally bound and the non-marginally bound solutions of 
the 5D perfect fluid collapse case [10]. 
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IV. APPARENT HORIZONS 



In this section, we discuss the apparent horizons and their physical signifi- 
cance by using the solution to the field equations. The apparent horizons can 
be found by using the boundary of the three trapped spheres whose outward 
normals are null. For the interior metric, this is given as follows: 

g^Yv = Y 2 — (^) 2 = 0. (60) 

Inserting Eqs. f l4"0]) and (jHJ) in the above equation, it follows that 

(A + K 2 - 87rp c )Y 4 - 6Y 2 + 12m = 0. (61) 

In particular, when we take A = 87rp c — K 2 , Y = \j2m. This is called the 
Schwarzschild horizon. For m = 0, po = 0, and K = 0, we have Y = \J\, 
which is called the de-Sitter horizon. The following positive roots are found 
from Eq. ( 16"TT) . 

Case (i): For 4m < A _ 87r p o+y2 , we obtain two horizons: 



n = \ - x — - — ^+ v v Z (62) 





a/9 


- 12m(A - 


8tv Po + K 2 ) 


A-87rpo + ^ 2 ' 




A - 8np 


+ K 2 






- 12m(A - 


87rp + K 2 ) 


A - 8np + K 2 




A - 87rp 


+ K 2 



Y 2 = \ — -- v ^~^ ^I^Ll . ( 6 3) 



When m = 0, these reduce to Y\ = y (a + k^-8tt P ) anc ^ ^ 2 = 9- ^i an d Y 2 
are called the cosmological horizon and the black hole horizon, respectively. 
For rn^O and A ^ 8np c — K 2 , Y\ and Y 2 can be generalized [20], respectively. 

Case (ii): For 4m = , 3 , we have repeated roots: 

V ' y/(A+K 2 -87rp c ) 1 

Y 1= Y 2 = 3 = Y, (64) 

x/(A + K 2 - 8np ) 

which shows that both horizons coincide. The ranges for the cosmological 
and the black hole horizons are 



° 5 * * ^ U-J n + K> £ * £ V A-8^ + JP - (65) 
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The black hole horizon has its largest proper area AnY = ( A+A . 2 _ 87rp \ , and 
the cosmological horizon has its area between ov+^^87rp^) anc ^ (A+K^-8np ) ' 

Case (iii): For Am > , 3 , there are no positive roots; conse- 

quently, there are no apparent horizons. 

Now we find the formation time for the apparent horizon with the help 
of Eqs. f )53|) and flBTj) given by 



t "^°-'/ 2(A-84.H-^) Smlr ' ( lr- 1)i - ( " =1 ' 2) - (66) 

In the limit (8717)0 — K 2 ) — > A, we obtain the result corresponding to the 5D 
Tolman-Bondi solution [19], 

tah = t - M. (67) 

Equations flB^]) and (jSHD imply that Y\ > Y 2 and t 2 > ti, respectively. The 
inequality t 2 > ti indicates that the cosmological horizon forms earlier than 
the black hole horizon. 

V. OUTLOOK 

This paper provides an extension of our previous analysis for an electromag- 
netic field in perfect fluid gravitational collapse with a cosmological constant 
[14] from the AD to the 5D case. The exact solution for the interior space- 
time with a charged perfect fluid in the presence of a positive cosmological 
constant is derived. The effects of an electromagnetic field on the 5D gravi- 
tational collapse follows: 

The relation for the Newtonian potential is 4> = |(1 — goo) . Using Eqs. 
f fT2|) and fj49|) . for the exterior spacetime, the Newtonian potential turns out 
to be 

0( J R) = ^ + (A-8 7 rp o + K 2 )^. (68) 
The corresponding Newtonian force is 

F = -^ + (A-8vrp + K 2 )|. (69) 
This force will be zero for 
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R = r and m = r while it becomes repulsive 

(A-8Trp +K 2 )? 12(A-8tt Po +K 2 )? 

(attractive) for larger (smaller) values of m and R. For the repulsive case, 
one must have A > (87rpo — K 2 ) such that 87rpo > K 2 over the entire range of 
the collapsing sphere. Notice that K = K(r) gives the electromagnetic field 
contribution. From Eq. (I44p . the rate of collapse turns out be 

2m Y 
Y = -— + (A-8np + K 2 )-. (70) 

Here, we have re-formulated the Newtonian model in terms of the acceleration 
of the collapsing process. If 8np — K 2 > A over the entire range of the 
collapsing sphere, then the force becomes attractive, and the cosmological 
constant would favor the collapsing process. 

Further, we have found two apparent horizons (cosmological and black 
hole horizons) , whose areas are larger due to the extra dimension as compared 
to the values obtained from our AD analysis [14]. The solution for the non- 
marginally bound case (F(r) ^ 1) is also possible in 5D, which is not possible 
in the AD case. Equation ( HHj) represents the total amount of the collapsing 
mass contained in a sphere of radius to r. This amount of mass is larger 
than that to the AD case; hence, collapsing process is faster in this case. 

We would like to mention here that all our results match to 5D perfect 
fluid case [10] if we take K = 0. However, the results for a AD gravitational 
collapse with an electromagnetic field [14] cannot be recovered directly. The 
reason is that Eq. (61) is quartic polynomial while in the AD case, the 
corresponding equation is a cubic polynomial. In the 5D case, the roots 
under the possible conditions are different than the roots in AD. Thus, there 
is only a difference of one degree in both cases, but the solutions are much 
different, depending on the nature of the polynomials. 

The cosmological horizon was found to form earlier than the black hole 
horizon. Also, Eq. ( 1661) shows that the apparent horizon forms earlier than 
the singularity; hence, the end state of the gravitational collapse is a black 
hole. These evidences agree with the AD case. There was a possibility of 
a locally naked singularity in the AD case due to the presence of an elec- 
tromagnetic field, but such a possibility may be avoided due to the extra 
dimension. Thus, we conclude that 5D favors the formation of a black hole. 
It is interesting to mention here that our study supports the CCH. 
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